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Abstract 

We derive the loop equations for the one Hermitian matrix model 
in any dimension. These are a consequence of the Schwinger-Dyson 
equations of the model. Moreover we show that in leading order of 
large N the loop equations form a closed set. 
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1 Introduction 



Recently enormous progress has been made in the understanding of the 
physics of non-perturbative strings in a low number of dimensions. This has 
been made possible by the matrix model technology. The Feynman graphs 
of the matrix models provide a suitable discretization of the two-dimensional 
surface where the string lives, classified according to its genus, which in the 
matrix model is a given order in 1/N (N being the range of the matrix). By 
tuning the couplings in the matrix model in such a way that the perturba- 
tive series diverges, the continuum limit of the two-dimensional surface can 
be approached. In this way it is possible to get a sum over all topologies 
in the first quantized string, using as a discretization of the two-dimensional 
world-sheet a dense sum of Feynman graphs of an associated matrix model. 

Whenever the matrix model is solvable, the corresponding string prop- 
erties can be also determined. In this way 2d gravity coupled to conformal 
matter with c < 1 has been analysed and the critical exponents determined 

I@ I- 

Furthermore it has been shown that a Virasoro algebra naturally arises in 
these matrix models This algebra is a realization of the Schwinger-Dyson 
equations in the U(N) invariant set of operators of the matrix model. Using 
this algebra, the equivalence of zero- dimensional matrix models to topological 
2d gravity has been proved [||. Some time ago, an application of these ideas 
to compute intersection indices on moduli spaces of Riemann surfaces has 
been made ||. 

More recently further progress has been made possible by the complete 
solution of the two-matrix model in zero dimensions 0] [§ ||] . 

To study non-critical strings in higher dimensions, a knowledge of the 
solution of matrix models in d > 1 is needed, d is the dimension of space- 
time, (for a different approach to the problem see |ITJ). But beyond d = 1 the 
matrix model is not solvable (although the gauged matrix model is [lW). So 
various approximation schemes have been suggested to tackle this problem 



In this note we want to derive the loop equations for the Hermitian ma- 
trix model for arbitrary d. They are a consequence of the Schwinger-Dyson 
equations that its correlation functions satisfy. Moreover we will show that in 
the large N- limit, the loop equations form a closed set. We believe this is an 
important step in the understanding of the physics of these models because 
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recent experience has shown [ffl that the Schwinger-Dyson equations are a 
powerful tool to explore the involved symmetries . It would be desirable to 
explore the loop equations we present here from this perspective. 

A simple and powerful way to derive Schwinger-Dyson equations is the 
hidden BRST method [pL3[| . In fact this was used to find the loop equations 
for the zero-dimensional, two- matrix model in [Q. So we start by reviewing 
this method in section 2 and apply it to derive the loop equations of the 
(i-dimensional matrix model in the next section. In section 4 we consider the 
case d = as a check. In section 5 we draw some conclusions and comment 
on future work. 



2 Hidden BRST Symmetry and Schwinger- 
Dyson Equations 

The common procedure to derive Schwinger-Dyson equations (SDe) involves 
using the invariance of the path integral measure under field translations. 
Invariances of the action itself are not relevant for this derivation. In this 
section we review a different approach which has been proposed some years 
ago [13| : By making use of a (trivial) BRST symmetric extension of any 
action S, we can derive all SDe's as BRST supersymmetric Ward identities 
of the new action. 

For the sake of simplicity consider the path integral describing one bosonic 
field <p(x): 

Z = j[dtp]e- s ^. (1) 

The functional measure [dip] in the last equation is invariant under field 
translations: 

<p(x)^<p(x) + e(x). (2) 
The invariance of the functional measure implies the following statement: 
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From (3) we get 

. SF SS . 

for any F = F((p). These are the SDe's for the theory. Notice that, in 
general, S is not invariant under (2) and that this lack of invariance plays 
no role in the derivation of the SDe's We now introduce auxiliary fermionic 
variables ip(x), iJj(x) and insert a trivial factor of unity 

[#][#]e - i 

inside the partition function. The resulting action S is invariant under the 
following transformation: 

S[<p,il>,$\ = S[<p] + J dxi>(x)i>(x) (5) 

5ip(x) = ip(x) (6) 

5ip(x) = (7) 

6$(x) = -8S/5ip(x) (8) 

Associated with the BRST-like symmetry, (6)-(8) is a set of Ward identities; 
for example, unbroken BRST implies: 

(5[F(<p)j(y)]) = 0, 

which is 



SF S S 

dx-—Mx)^{y) - F —) = 0. 
6(f(x) dtp 



Computing the average with respect to the fermionic variables we recover eq. 
(4). Notice that the BRST transformation we have just defined commutes 
with any symmetry the original action S may have. In particular, if the 
original action is invariant under a U(N) group, the BRST symmetry will 
commute with this group also and then its implications will be true order by 
order in the 1/N expansion. 

In the next section we use the Schwinger-Dyson BRST symmetry to derive 
the loop equations for the ci-dimensional Hermitian matrix model. 
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3 Loop Equations 

In this section we will derive the loop equations for the model defined by the 
action: 

S = J d d xtr (^d lx M(x)d^M(x) + V(M(x))j (9) 

M(x) is an N x N Hermitian matrix defined on the Euclidean space-time 
point x, and V is a local function of M (normally a polynomial whose degree 
can be restricted by renormalizability of the model for d > 2). 
The extended action and the BRST transformation are: 

S ext = S + J d d xtr$(x)^(x) (10) 
5M(x) = i/j(x) (11) 

5ip(x) = (12) 

WW = -jik ( 13 ) 

Introduce the following representation of M(x): 

N 

M(x) = Y,™ a (x)T a (x), (14) 

a=l 

where the projectors of the matrix M(x) satisfy: 

N 

J2T a (x) =1 (15) 

Q = l 

trT a (x)T^x) = 5 a(3 (16) 
T a (x)T^x) =6 aP T a (x). (17) 

If S denotes the BRST variation (see last section), we get the following BRST 
transformation for the projectors T a and eigenvalues m a (see the Appendix): 

ST a (x) = E/Ma ^W^jxj+^W^Cx) (lg) 
v / i—ip^po. m a (x)—mp(x) \ ' 

5m a {x) = trip{x)T a {x). (19) 

Consider the following functional of X(x) 

N 

K tj = dx^K(x)T a (x) ir (20) 

J a=l 
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Notice that Kij is an ^-independent N x N matrix. 

The basic object we will use to write the loop equations is: 

u[X\ = e lK (21) 

It fullfils the identity: 

8u[X] 

5Xjx) Jo 

which implies: 



/ dte ltK iT a {x)e t{1 - t)K , (22) 

JO 



The trace is taken with respect to the internal indices only. From u[X] we 
can compute all symmetric combinations of products of T a . For instance: 



I A=0 



The Schwinger-Dyson equation is: 

{6tr[u[\]$(y)]) = (25) 

i.e. 

(tr[5u[X]ij(y) + tr[u[\]6ip(y)]) = 0. (26) 

Since 

6$(y) = U y M - A(M) (27) 

where 

A(M) = V'(M). (28) 

We get: 

^[□ 2/ M(y)-A(y)] = J2l m M D y + 2 9 ll m a (y)d ll + a y m a (y) - A a ]truT a (y). 

a 

(29) 

We can evaluate: 



We can also compute 5u [A]: 



jQ •* a=l p^a 

T a (x)iP(x)T?(x) +TP(x)^(x)T a (x) cl(1 _ t) j dx ^ K{x)Ta{x) 
m a — mp 

After computing the fermionic average, we get: 

<<W(V)> = (i £ t K{V) ~ tru[t\}T\y)tru[{l - t)X]T a (y)). (32) 
p^Jo m a — mp 

Since T a (y) appears only once in the trace, we can express it as a deriva- 
tive of u: 

M(1 _ t)A]n ,)^^=M (33) 
Therefore the loop equation is: 

ftfaJo m a -m p t{l - t) 5Xp(y) SX a (y) 
-(E[m a (yPy + 2d^m a {y)d^ + U y m a {y) - A a ]^^), (34) 

which is valid to any order in the 1/N expansion. 

A simplification is possible if we restrict ourselves to the leading order 
in 1/N . In this case we can apply Witten's factorization property |T4]|to the 
loop equations to prove that: 

F(M) ~ (F(M)) + 0(N~ a ), a>0 

for any [/(iV)-invariant F. In particular, we obtain that: 

m a (x) ~ (m a (x)) = (m a (0)). 
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Therefore, in leading order of large N the loop equation becomes a closed 
system for u[X]: 



K{y) ~ Xp(y) 1 Stru[t\] Stru[(l - t)X] 



^ a Jo m a -mp t(l - 1) 5Xp(y) SX a (y) 

"E^-^Hf < 35 > 

The reason why we could derive this closed system was that we chose to 
write the Schwinger-Dyson (loop) equations in terms of T a rather than M. 
Notice also that by using the BRST Schwinger-Dyson transformation we can 
avoid this difficult change of variables in the path integral and the (highly) 
non-trivial measure in the (m a ,T a ) basis. 

The loop equations must be solved subject to the following boundary 
conditions: 

tru\ x =o = N 
Stru 



i 



5X a (x) 

^ tru i x aa\ 

r X / x x / x U=o - -o a p {J® 
5X a (x)Xp(x) 

A particular solution of the loop equations is u — constant (A indepen- 
dent), which does not satisfy the boundary conditions. We must look for a 
non-trivial solution. 

An equivalent formulation of the loop equations is obtained by noticing 
that only the derivative of tru appears in the loop equations. So identify: 

v[X]a(x) = TT~7~T- (37) 

oX a {x) 

Then the loop equation is equivalent to the following system of equations: 

E I' dt K t ] ~ t {V) tn 1 r A txUy)v[(i-t)x] a (y) 

p +a h m a -mp t(l-t) 

= -J2(m a a y -A a ))v[X} a (y) (38) 
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SX (y) " 5X a (x) { ) 

v[X] a (x)\ x =Q = i (40) 
Sv[X] a {x) l 

-da/3 (41) 



5Xp(x) A 
Notice that: 

trM{x 1 )M{x 2 )...M{x n ) = ]T m ai ...m a JrT^{x 1 )...T a -{x n ) 

ai...a n 

= — m ai ...m an [trT ai (xi)..T an (x n ) + all permutations of indices a] (42) 

ai...On 

Therefore it is enough to know trw[A] to get all correlations of the matrix 
M that are symmetric under permutations of the space-time points X{. 

We can get a continuum version of the large- N loop equation by the usual 
manipulations. That is we introduce the functions X(x, y) and m(x) defined 
for < x < 1 and the density of eigenvalues p, and make the following 
identifications: 

X a (y) = X(^,y) (43) 
— ct 

m a = ViVm( — ) (44) 

Km) = ^ (45) 
p(m) = 1. (46) 

Then the large- iV loop equation becomes: 

^ l X(m 1 ,y)-X(m 2 ,y) 1 



// 
dm 1 dm2p(m 1 ) p(m 2 ) J dt- 



m\ — m 2 t(l — t) 
, 5tru[tX] 8tru[(l - t)X] = _ r _ tru[X] 

5X(m 2 ,y) SX(mi,y) J 5X(m,y) 



P stands for the principal value of the integral over m,. 

In the next section we will solve in detail the d = case and show that 



it agrees with the known solution [15 
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4 Zero-dimensional matrix model 



The zero- dimensional matrix model was solved a long time ago [15]. In this 



case the only unknown variable is the density of eigenvalues p. This model 



was solved using the SDe's in [16 



It is easy to show, using the properties of the projectors, that: 

u[\] = Y,e iXa T a (48) 

a 

tru[X] = J2e iXa (49) 

a 

That is, in d — 0, trit[A] decouples from M. Replacing trit[A] into the loop 
equation we get the equation that determines the eigenvalues: 

[2£ A J =0 (50) 

Since X a is arbitrary, we obtain the known answer: 

A a = 2E " (51) 

5 Discussion and Open Problems 

We have been able to derive the loop equations for the Hermitian matrix 
model in arbitrary space-time dimension d in terms of a functional it [A]. 
Moreover, in leading order of large N, the loop equations form a closed 
set. All this was possible because we had chosen to write the SDe's in the 
projector-eigenvalue basis instead of the M-basis, as usually done in zero 
dimensions. 

We believe this to be an important step in the comprehension of the 
physics of this model. Although no solution of the equations is available at 
present, a numerical solution may be feasible. Also a variational ansatz of 
the type employed in |T7| for the QCD string may work. 



From a conceptual point of view, our loop equations may shed light into 
the geometrical structure of the model. The <i-dimensional matrix model 
provides a discretization for a string in d + 1 dimensions, but the exact 
geometrical mapping between the two descriptions has not been done yet. 
As in QCD |18[ may be the loop equations will help to bridge such a gap. 
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Appendix 

In this Appendix we want to derive the Schwinger-Dyson BRST variation 
of the eigenvalues and projectors of the matrix M. Since the space-time 
dependence of the involved variables is irrelevant to this discussion, we will 
not write it explicitly. 
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To begin with, we write the eigenvalue equation satisfied by the projec- 
tors: 

MT a = m a T a . (A.l) 

Computing the BRST variation of this identity, multiplying by T 13 from the 
left and symmetrizing the resulting expression in a and (3 we get: 

T a i)T p + T^T a = (m a - mp) [T p 5T a - T a 5T p ) + 25m a T a 5 a p. (A.2) 

From a = f3 we get 

T a ipT a = 5m a T a . (A.3) 

That is 5m a is the eigenvalue of T a ip of eigenvector T a . Computing the trace 
of the last equation we get: 

5m a = tripT a . (A.4) 

Consider now the case a 7^ (3 in (|A.2 ). Dividing by m a — mp and summing 
over (3 7^ a we get the BRST variation of the projectors: 

«- = v I!^±I!^!. (A . 5) 

As a check, we can easily show that these transformation rules for the 
projectors and eigenvalues of M(x) are compatible with eqs. (15)-(17) and 
that they reproduce eq. (11) as well. 
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